The symbol of a function of an operator 

Alfonso Gracia-Saz * 
February 1, 2008 

Abstract 

We give an explicit formula for the symbol of a function of an opera- 
tor. Given an operator A on L 2 (R N ) with symbol A G C°°(T*R N ) and a 
smooth function /, we obtain the symbol of f(A) in terms of A. As an ap- 
plication, Bohr-Sommerfeld quantization rules are explicitely calculated 
at order 4 in h. 



1 Introduction. 

The goal of this paper is to give a realistically computable formula for the 
symbol of a function of an operator. Let A be an essentially self- adjoint operator 
in L 2 (R N ) with symbol (definitions to follow; see A e C°°(T*R N ). Let 
/ : R — ► R be an analytic function. Let B = f(A) be an operator with symbol 
B. We want to write B in terms of A. We will derive the formula: 

The sum is taken over finite graphs Y with no isolated vertices. For every such 
graph r, V is the number of vertices and E is the number of edges. Ar(^4) is a 
polynomial in the derivatives of A constructed algorithmically from T (see © ■ 
Sr is the order of the symmetry group of T. cr is a simple invariant of T (see 
H4.3J1 . The terms through order 4 in ft of l|l.lf) are shown in Appendix iDl 

The existence of a universal equation like Hl.lfl was used by Voros an( A 
Colin de Verdiere [H] as part of a calculation to obtain Bohr-Sommerfeld quan- 
tization rules at higher orders in H. They derived it by using a explicit spectral 
theorem that writes f(A) in terms of the resolvent. 1 Their method gives a recur- 
sive way of obtaining higher order corrections in ft but is, in practice, intractable 
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1 In Appendix [E] we give an alternative derivation of our main result inspired by this 
approach. 
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after order 2. In contrast, the diagrammatic notation that we use (inspired by 
[J) makes it simple to derive all our formulas and to write down explicitly their 
terms. 

For the derivation of (|l.lf> we need: 

• Weyl quantization, i.e., a well defined correspondence between operators 
and symbols (see M2.1|l . but we will not use its explicit form. 

• The explicit form of the Moyal product (sec H2.2I) . 

• A spectral theorem, i.e., a way to define a function of an operator (see 
^2.3fl . but we will not use its explicit form. 

Before going on, let us mention some possible applications of this calculation: 

1. Bohr-Sommerfeld quantization rules. This is treated in |Q 

2. Determinant of certain differential operators. See |14| . Such determi- 
nants naturally arise in quantum field theory at the one loop level. The 
determinant of an operator can be defined with the property log det A — 
trace log A, which holds in finite dimensional spaces. Then, the trace of 
the operator B — log A can be calculated by means of its symbol B. 

3. Star exponential of quadratic forms. Omori et al analyzed this problem in 
|16| . In the case when A is a quadratic function and / is an exponential, 
formula simplifies. See 313 

The structure of this paper is as follows. In &\ we explain the necessary 
background: Weyl quantization, the Moyal product and functions of operators. 
In |21 we introduce diagrammatic notation that will be used afterwards. ^1 i s 
the core of this paper, where we derive and other equivalent equations. 

In fJSJwe consider generalizations to functions of various variables and to other 
quantizations. In SjSJwc study the case of a quadratic symbol, in particular the 
harmonic oscillator hamiltonian. This is an example of restricting to a smaller 
class of symbols A, for which the family of graphs to consider becomes smaller, 
too, and simplifies. |J7|explains the application to Bohr-Sommcrfeld rules. 
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thankful to Matthew Cargo, Laurent Charles, Yves Colin de Verdiere, Robert 
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2 Background. 



2.1 Weyl quantization. 

Weyl quantization [201 (or Weyl-Wigner correspondence) is a way to relate the 
classical and quantum descriptions of a system. In the classical description of a 
system, the space of states is a Poisson manifold, whereas the quantum space is 
a Hilbert space. For a particle in N dimensions, the Poisson manifold (classical) 
is T*R N and the Hilbert space (quantum) is L 2 (R N ). 

The observables are classically described by smooth functions on the phase 
space C°°(T*R N ). In the quantum description they are operators on the Hilbert 
space Op(L 2 (R JV )). For instance, we have canonical coordinates 

(x,p) := (xi, . . .,x N) pi, . . . ,p N ) 

on the cotangent bundle T*R N : Xj are the coordinates on the base R N and 
Pj are the coordinates along the fibers. Xj and pj are elements of C°°(T*R N ) 
and hence classical observables. We associate to them operators Xj and fj on 
L 2 {R N ) defined by: 

Xj[<p](x) —Xj(p(x) 

fj[4>](x) = - (x) 

for G L 2 (R N ). 

Here H is the Planck constant. We will treat it as a formal parameter and 
consider that all our spaces are formal power series in h. Thus the classical 
observables will be elements of C 00 (T*R Ar )[[?i]] and the quantum observables 
will be elements of Op(L 2 (R JV ))[[ft]]. 

A quantization is an extension of this correspondence to a map 

AeC°°(T*R N )[[h}} ^ AeOp{L 2 (R N ))[[h}} 
In particular the Weyl-Wigner correspondence is defined by: 

In order for it to be a well-defined bijection we need to restrict the domain 
and consider only a certain family of smooth functions and a certain family of 
operators. See 0, for details. In H5.2l we consider alternative quantizations. 

A is called the symbol of A. Moreover, the space of operators with compo- 
sition is an algebra. We can define an associative operation in C°°(T*K A ')[[fi,]], 
called star product *, that makes the bijection A i— > A into an algebra isomor- 
phism. In other words, C * D is the symbol of the operator CD. 
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2.2 The Moyal product. 

Moyal [TS] gave an explicit expression for the star product in the case of the 
Wcyl quantization, called the Moyal product (but actually due to Groenewold 
|3 )■ It is derived from the definition of the Moyal product (C * D = CD) and 
the explicit form of Weyl quantization |[OJ|. If C, D are symbols in C°°(T*R N ), 
then: 

We need to explain the notation in 12.2|l . In the natural coordinates (x,p) on 
T*R N , the Poisson bivector field is 

N 

J = ^2(d Xj ® d Pj - d Pj ® d Xj ) 

3 = 1 

where d q := J^. Let us call the coordinates (z 1 , . . . , z 2N ) := (x,p) and J^ v the 
coefficients of the Poisson tensor on this chart 

where In is the identity N x N matrix. Using Einstein summation criterion 
(summation over repeated indexes), J = J^d^ <g> d v . Let us also write 

dC 

C, M = d^C = d z »C = — 

d k C 

C^i.-./ij. = . . . d fJik C = d z n . . . d z v k C 



dz^ . . . dz^k 

Then the Poisson bracket in T*M. N can be written as 

{C,D} = C t ^C !V 
The terms appearing in l|2.2[) are defined as: 
{C,D} = CD 

{C,D} 1 = {C,D} = C. p J fW D^ 

{C, D} k — C^i . ^fc J M v ■ ■ ■ J M D v i v k 

Following 0] we use the notation " — >" for a Poisson tensor in the following 
way: " — >" is replaced by J M1/ , the expression in the head of the arrow is acted 
on by d v and the expression in the tail of the arrow is acted on by d^. For 
instance 



C =C^J^D, V = {C,D} 

C ?D =C^^J^^J^- 2 D^2 = {C,D} 2 
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Here, \k denotes k arrows. The same can be done with more complicated dia- 
grams: 



JM2"2 

C — D ^ A = C^D^A (23) 

u ,/il ,y ^,^1/^2/^3 ^ L ,^2I-'3 

Since J MI/ is skew-symmetric, inverting an arrow multiplies the expression 
by -1: 

C — >D = (-1)C< — D 
With this notation, the Moyal product (|2.2|l is written: 

The fact that the Poisson bracket is a derivation on each argument (Leibniz 
rule or product rule) is written as: 

{CD) -,E = C {D^E) + (C^E)D 
{CD,E} = C{D,E} + {C,E}D 

Another application of the product rule is 

(C ->■ D) E C ->D ->E + E^C D 

(C^J^D tl/ ) >a J^E <0 = C.^D^J^E.p + G^J^D^E^ 

(2.5) 

Notice that (|2.5|l cannot be written with Poisson brackets { , }. 

This notation makes it simpler to write certain calculations. Lemma I A. II in 
Appendix lAl is a generalization of (|2.5() that will be used in our derivations. 



2.3 Spectral theorems. 

To understand we need to define what a function of an operator means. 
If f(y) — y n then f(A) — (A) n . A spectral theorem extends this definition 
of f(A) to a wider class of functions /. More specifically, let A be an algebra 
of smooth functions under pointwise multiplication and let A be an operator. 
Then a spectral theorem is an morphism of algebras 

/eiH f(A) e o P (L 2 (R N )) 

with certain properties. See \7\ for details. There are many spectral theorems 
(i.e. different algorithms to calculate f{A)) for different algebras A. See pQ for 
a complete list of references. We will not use any explicit form of a spectral 
theorem (with the exception of Appendix [51 where we derive an alternative 
proof of our main result). We just need the fact that a spectral theorem is a 
morphism of algebras: 

f(A)g(A) = (fg)(A) 
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3 Graphs. 



The formulas we are going to derive are power series whose terms are labeled 
by graphs. We define now the family of graphs that we are going to use and 
introduce notation. 

A graph consists of a finite set of vertices and a finite number of edges. Each 
vertex is represented by a dot. Each edge is represented by a line joining two 
vertices. Multiple edges joining the same pair of vertices are allowed. A self- 
edge (an edge from a vertex to itself) is not allowed. A graph does not need to 
be connected. 

An example of a graph is 

• » • = • • • (3.1) 

A labeled graph is a graph with in which we have labeled the vertices with 
the first V natural numbers 1,2, ... ,V and the edges with the first E natural 
numbers 1, 2, . . . , E. Therefore, V is the number of vertices and E is the number 
of edges. For instance: 

2 

1^^2^^6=f=3 4^^5 

Even though a graph as defined above is not oriented, a labeled graph has a 
natural orientation: every edge is oriented so that the target has a higher label 
than the source: 

2 

1^^2^*6*1=3 4^-5 

All the information in a graph is given by how many edges there are joining 
each pair of vertices. All the information in a labeled graph is given by which 
edge joins which pair of vertices. This leads us to adopt the following formal 
definitions. 

Definition 3.1. Let V and E be two non-negative integers. A labeled graph 
with V vertices and E edges is a map 

s:{l,...,E}^V 2 {l,...,V} 

where ViX denotes the set of subsets of X with 2 elements. 

This means simply that the edge i joins the pair of vertices s(i). The group 
Se of permutations of E letters acts on {1, . . . , E}. The group Sy of permu- 
tations of V letters acts naturally on P%{1, . . . ,V}. Hence the direct product 
Sy x Se acts on functions {1, . . . , E} — ► Vi{\, ■ ■ ■ , V}, that is, on the set of 
labeled graphs with V vertices and E edges. 

Definition 3.2. An unlabeled graph or simply graph with V vertices and E 
edges is an orbit of this action. 



6 



We will denote a labeled graph by T and the corresponding (unlabeled) 
graph by [r], if wc need to distinguish between them. Otherwise, we will abuse 
notation and denote a graph simply by T. 

It is convenient to define now two more concepts that will be needed later. 
The order of symmetry of a labeled graph is the number of permutations of 
edges and vertices that we can make without changing it. Or, more formally: 

Definition 3.3. The order of symmetry St of a labeled graph T is the order of 
the stabilizer of T in the action of Sy X Se on the set of graphs with V vertices 
and E edges. 

See Appendix IU1 for examples. 

A labeled graph is reduced if it does not have any isolated vertices. Or, more 
formally: 

Definition 3.4. A labeled graph s : {1, . . . , E} — ► T^lL ■ ■ • , V} is reduced if 
every i = 1, . . .V is in some element of the image of s (i.e., "if every vertex is 
in some edge"). 

Both concepts (order of symmetry of a graph and reduced graph) extend 
naturally to unlabeled graphs. All the previous examples are reduced. The 
graph • • • is not reduced. 

Given a labeled graph T with V vertices, and given V symbols Ax,..., Ay 
we construct a new symbol, called Ar(^4i, ■ • • , Ay), by substituting the vertices 
with Ax,. . . Ay, and letting every edge represent a Poisson tensor (as explained 
in ESI- We denote X r (A, . . . , A) simply by X r (A). 

For instance, if T is the labeled graph: 




then Ar(C, D, A) is the expression in 1|2.3[) . And 

X r (A) = A »■ A *- A 

= A »- A^=*A 

JM3"3 

_ A TMlfl A A 

— ■"■,111" rl ,I'lA'2M3 ° ° ■ rL ,V2V3 

Since changing the direction of one arrow multiplies the expression by — 1, 
Ap] is defined up to a sign. 

4 Main results and calculations. 

We recall our problem. Let us fix an operator A with symbol A and a smooth 
function /. Let B = f(A) be an operator with symbol B. In this section we will 
perform the necessary calculations to obtain various expressions of B in terms 
of A. 
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The main step is to obtain an expression for an iterated star product C\ * 
. . .-kC n for symbols Ci in terms of graphs. We do this in H4.ll Then in H4.2l we 
derive our hrst expression for B in terms of A. 

Equation l|4.4|l at the end of SJ4.2l is a power series whose terms are parametrized 
by labeled graphs. This is the easiest form of our result to derive, and it is use- 
ful for theoretical proofs. However, it is not convenient for explicit calculations 
when we want to write the first few terms explicitely. There are only a few unla- 
beled graphs, but many labeled graphs. In H4.3l we obtain our second expression 
for B in terms of A, 1)4. 8|) . a series whose terms are parametrized by unlabeled 
graphs. There is still a third form of our formula, Equation 14.12fl . whose terms 
are parametrized by connected graphs. This last form is studied in H4.4I 

Using either of these equations, we have included in Appendix ID1 the explicit 
form of the terms up to order 4 in h of the symbol B in terms of A. 

4.1 The n— th star product. 

The main step in the derivation of l|l.lfl is the following expression for the 
iterated star product, which generalizes Moyal's formula: 

Lemma 4.1. Let C\, . . . , C n 6 C°°{T*M. N ) be symbols. Then 



Note that when n — 2, the previous lemma is exactly the Moyal formula 
for the star product (|2.4II . For every k, there is only one labeled graph with 2 
vertices and k edges. 
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In the previous expression, the vertices are labeled from left to right 1,2,3 in 
all graphs. 

Proof of lemma \^~l\ We will use induction on n. 
The result is true for n = 0, 1, 2. 

Inductive step. We use the associativity of the star product. 
C\ * . . . * C n+ i = (Ci * . . . C n ) * C n+ i 



labeled 
graphs r 
with n vertices 



E E 

labeled fc— 
graphs r 
with n vertices 



n+1 



i an 



E\k\ \ 2 



E+k 



Ar(Ci, ■ ■ ■ ! C n ) \k C n +i 



We can apply lemma HOI in Appendix 1X1 to Ar(Ci, . . . , C n ) \k C„+i 
and we will get a sum over labeled graphs V with n + 1 vertices! Triey are 
built by starting with a labeled graph T with n vertices (labeled 1, . . . , n) and 
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E edges (labeled 1, . . . , E), adding the (n + l)-th vertex and k edges (labeled 
E + 1, . . . , E + k) ending at the (n + l)-th vertex. The number of edges of T' 
is E 1 — E + k. If we want to account for all possible labeled graphs with n+1 
vertices and E' edges, we need to divide by a factor of ( E ^ k ) in order to account 
for all ways of relabeling the edges. Fortunately, we can write 

1 1 fE + k\ 

EW ~ (E + k)\ V E ) 

and we get 

C\ ★ . . . * c„+i = 

labeled 
graphs r 
with n+1 vertices 

□ 




4.2 First formula for the symbol of a function of an oper- 
ator. 

We now attack the problem of obtaining the symbol B of B — f(A) in terms of 
A. 

To calculate B we use the following expression [5]: 
Lemma 4.2. Let z e T*R N and let a Q := A(z ). Then 

B ^ - E y\ f{k){ao) (A - ao) * k {zo) (4 - 2) 

fc=0 

Note that the right hand side of l|4.2|l is well-defined only at the point z . 
To prove it, we need the following fact: 

Lemma 4.3. Let g be a smooth function and C := g(A). Let oq := A(zq). If g 
has a zero of order m at oq, then C(zq) — 0(h m / 2 ). 2 

Proof of lemma \4-3\ Let us write g(y) — gi(y)(y — a ) m and let C\ = gi(A). 
Then C = g(A) = gx{A){A - a ) m and C = Ci* (A — a )* m . Since (A(z) - 
ao)\ z=zo = 0(h), then C(z ) = 0{h m ' 2 ). □ 



Proof of lemma \J~^ For every m, apply lemma 11751 with 

m 

g(y) :=f(y)-J2uJ ik) ( a °)(y- a °) m 

k=0 

This proves that 

B ^ - E h f(k){ao) ( A - a oT k ( z o) =o(n^ 2 ) 

k=Q 

for all m. □ 



2 Actually, we can do better: C(z ) = 0(?i 2m/3 ). 
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Now we only need to substitute into l|4.2|l to get an expression for B 

in terms of graphs: 

co co , / j- \ k 

= E n\ f(n){A{ZQ)) E E fci ( V) A r(^-«o)(.o) 

n— fc=0 labeled graphs T 

with n vertices 
and k edges 

(4.3) 

In order to calculate Ar(^4 — ao)(zo) we need to write the symbol A — ao at 
every vertex of T. If a vertex is not isolated, then some derivatives are acting 
upon that symbol, and we may substitute A— ao with A, since ao is a constant. If 
a vertex is isolated, then it contributes a factor of (A — ao), and (A — ao)(zo) = 0. 
Hence we only need to consider graphs without isolated vertices, which we called 
reduced. Putting this all together: 

reduced labeled 
graphs r 

It is to be noted that the previous equation is not, strictly speaking, a power 
series expansion in h, as A itself depends on h. However, it is very simple to 
write A = J^k h k Ak as a power series in h and expand (|4.4I) . As a matter of fact, 
we have an alternative way to write (|4.4|l with terms parametrized by graphs 
with weights where every term is a monomial in h. We believe that the approach 
shown in this paper is simpler, though, and clearly illustrates the method. 

Equation l|4.4(l was easy to derive, and it is useful for proofs and theoretical 
calculations, as well as to generalize to other quantizations (see tJ5.2[l . However, 
when we want to explicitly write the first few terms of it, this is not yet our ideal 
expression. We can put together the contribution of labeled graphs that differ 
only in the labels to get a series whose terms are parametrized by unlabeled 
graphs. We will do it next. 

See Appendix 151 for an alternative derivation of 1)4. 4JI . 



f [V> (A) 



VI 



Xr(A) 



(4.4) 



4.3 Version with non-labeled graphs. 

Two labeled graphs which are the same except for the labeling of vertices and 
edges give the same contribution (up to a sign). 
1)4.4(1 can be rewritten as: 

reduced v ' ' T'e[r] 

graphs [r] 

In words, we need to sum Ar" when T' runs through all possible relabelings of 
r. Remember that Ar' and Ar will be equal up to a sign. 
Define 

^ ^ ^ ^ 1)^ °^ arIOWS i nver t c d g°i n g from T to F ^ 

r is a reordering of 
the vertices of T 
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That is, we start with a labeled graph T. Then, we consider the VI possible 
ways of numbering the vertices of the graph with 1, 2, • • • , V. For each of them, 
we orient the arrows so that they all go from the vertex with the lowest label 
to the vertex with the highest label. Then we count these V\ relabelings with a 
sign, depending on the parity of the number of arrows inverted from our original 
orientation. 

We can then write the contribution from (|4.5|l as: 

J2 \ r <(A) = ^c r \ r (A) (4.7) 
r< e[ r] * 

where Sr is the order of the symmetry group of the (unlabeled) graph [r] (sec 
Definition 13. 3|) . The contribution corresponding to different relabelings of the 
edges is El. The contribution corresponding to different relabelings of the ver- 
tices is in cr- And we have to divide by the order of the symmetry group, to 
account for the situation in which exchanging edges or vertices results in the 
same labeled graph. 

For instance, if T is the graph then the contribution from 

renumbering the vertices is: 

1 ^ 2 ^ 3 1 ^ 3 ^ 2 2 ^ 1 ^ 3 

c r = (-1)° + (-1) 1 + (-1) 1 



2 -> 3 <- 1 



3 <- 1 -> 2 



3 <- 2 <- 1 
(-1) 2 



Finally we just have to substitute 14. 7|) into (|4.5|l : 



B = 



E 

reduced 
graphs [r] 



S r V\ 



\t{A) 



(4.8) 



Notice how cr and Ar(^4) are only defined up to a sign for [r]. However, 
those signs cancel in their product crAr(^4), which is well defined. 

The explicit calculation of Cr for a particular graph is actually very easy. See 
Appendix El which includes the value of cr and Sr for all reduced, connected 
graphs with 2 or 4 edges. Thanks to lemma IC.ll we only need to consider 
graphs r where every connected component has an even number of edges, since 
otherwise cr = 0. 



4.4 Version with connected graphs. 

Let us rewrite (|4.8|) as: 



B 



E 

reduced 
graphs [r] 



cr 
Sr 



Xt{A) 



VI 



f(A) 



(4.9) 
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where D is the differential operator which applies to /. 

Whenever we have an expression like l|4.9() . a series labeled by a certain 
family of diagrams, it is standard to reduce all calculations to only connected 
diagrams. Let's generalize. 

Let Q be the free commutative monoid generated by the set Go- Let S be 
a commutative ring (with multiplicative notation). Let O : Q — > S be a map 
satisfying 



0(r lXl 
for distinct x%, 



ni 



G N. Then, formally: 



E 0(x) = exp 



E °w 



(4.10) 



(4.11) 



In particular, consider Q to be a family of diagrams closed under topologi- 
cal sum and generated by the connected non-empty diagrams Go- If we write 
0(T) = , where Sr is the order of the symmetry group of the diagram and 
M. is a multiplicative function 

M(xi +X2) = M{xi)M{x2) for all xi,x 2 € G 

then O satisfies l|£TU|) . 
In our case, 



G = reduced graphs 

Go = reduced, connected, non-empty graphs 

5 , = C°°(T*IR JV ) 

■^([r])=(f) £ C r MA)^ 

M. is multiplicative from l|C.l|) in Appendix IO 



VI 



£ri frv 
V x ! "' K! 



Hence, using l|4.11|l in (|4.9|) : 



B 



exp 



E 



connected, reduced, 
non-empty graphs [r] 



V 



(4.12) 



5 Generalizations. 



As we mentioned in the introduction, we needed three things to derive our 
results: 



13 



• Weyl quantization. 

• The Moyal product. 

• A spectral theorem. 

Of these, we only used the explicit form of the Moyal product. Actually, the 
form of the Moyal product is calculated from the form of the Weyl quantization. 
We discuss now how to generalize to functions of various variables and other 
quantizations. 



5.1 Functions of several variables. 



Let A\, . . . , A n be n commuting operators in L 2 (R N ) with symbols A±, . . . , A n . 
Let F : M™ — > M be a smooth function. We consider the operator B = 
F(Ai, . . . ,A n ) with symbol B. Can we extend our results to calculate B in 
terms of A\, . . . , A n l The answer is yes. We are still using Weyl quantization 
and the Moyal product, and spectral theorems behave equally well for functions 
with several variables. See 1 for a list of references. 

Hence, we only need to repeat our calculations, starting from M4.2I but with 
a function of a several variables. The counterpart of (|4.8|l is 



B 



E 



reduced 
graphs [r] 



£r di % 
St 



d lv F{A) 



v\ 



(5.1) 



For instance the first few terms of (|4.8|l are 



B = f(A) 



4 



1 



.4: 



A 



f"(A) 



2 ^ 2! 
and the first few terms of l|5.1|l are 
B 



A^A^A 



f'"{A) 
3! 



F(A) 
4 



1 



A, 



.A, 



djdjF(A) 
2! 



A-, ► A* 



A, 



didjdkFiA) 



3! 



+ 0(ft 4 ) 

There are also the obvious versions with labeled or connected graphs. 



5.2 Other quantizations. 

There are other quantizations apart from Weyl quantization, that is, correspon- 
dances between operators and symbols |llj . If we want to use them, then we 
have a different star product instead of the Moyal product. The explicit form 
of the Moyal product has been used in two places: to prove lemma FOl and to 
derive an expression for the iterated star product (lemma |4.1|) . 



14 



Let us consider a generic star product that has the form 

oo -. / \ k 

C * D = E^( l i) {C,Dh (5.2) 

As long as (C, D) \— > {C, D}k is a bidifferential operator of order rrik and 
limfe_, 00 mj. = oo, then lemma 14.31 is satisfied. Hence, for those star prod- 
ucts, we only need to obtain an analogue to lemma FTTI that is, an expression 
for the iterated star product G\ * . . . ★ C n in terms of diagrams. This can often 
be done by induction if we start by writing the star product of two symbols as 
a series in terms of diagrams. 

For instance, we can consider standard order quantization: 

A(x,p) = u-p(x)A[u p ](x) 

where u p (x) = e lp ' x / h . In that case, the star product has the form of 15. 2|) with 
{C, D}k := dp.Cd^jD. Lemma T4. 31 still holds and, if we change the meaning of 
" — ►", lemma |4~T1 is also true. The simplest way is to use the same definitions 
we gave in H2.2I and Sjllto construct a polynomial Ar(^4) from a labeled graph 
r, but using the (non-Poisson) tensor 

(jn = ( I ° N °o) 

With that convention, (|4.4|l is still valid. However, it is no longer true that 
C — ► D = — C < — D. As a consequence, the versions with unlabeled or 
connected graphs are messier (although they still exist). 

In |13j , Kontsevich gave a star product that quantizes any Poisson structure 
on K . His expression is already a power series in h whose terms are labeled by 
a family of diagrams. We can use it to obtain an equivalent of lemma FTTI and 
we can derive, again, the counterparts of Equations Q4.4[l . I|4.8|l and l|4.12|l . In 
fact, our set of labeled graphs is a subset of Kontsevich's set of labeled graphs. 
It is to be noted, though, that Kontsevich's star product includes a weight 
associated to every diagram T which is, in practice, hard to calculate. (See [T2] 
and J7| for some results.) 

6 The case of a quadratic symbol. 

When we restrict to a smaller class of symbols, it is possible that the contribution 
of many graphs vanishes, simplifying our calculations. As an example, we study 
here quadratic symbols. 

A particular (simple) case of importance consists of taking A = I, the 
harmonic oscillator hamiltonian in 1 dimension. In the standard coordinates 
(z 1 ,^ 2 ) = (x,p) its symbol is I = ^((z 1 ) 2 + {z 2 ) 2 ). This simplifies the calcula- 
tions because any third derivative vanishes: I,^^^ = 0. The same is true for 
any quadratic function. Let be a 2 x 2 real symmetric matrix and consider 
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the symbol A = hz^QauZ" '■ Assume A is the symbol of an operator A. Then: 



z Q ^,1/ 

and any third derivative vanishes. As a consequence, we only need to consider 
graphs where every vertex has at most two edges. After lemma IC. II in Appendix 
1U1 we only need to consider graphs with an even number of edges. If we also 
ask them to be connected, reduced and non-empty, there are only two families 
of such graphs: 



A k = 
A fe = 



(2k edges) 
(2k edges) 



k > 1 
k > 1 



In words, A k consists of 2k + 1 vertices and 2k edges joined forming a line. 
Its symmetry group has order 2. A^ consists of 2k vertices and 2k edges joined 
forming a simple cycle. Its symmetry group has order 4fc. 

The corresponding polynomials in the derivatives of A can be calculated: 



— A 77*1 fl A TP.2V2 

= z^((QJ) 2k Q)^ 2k z^ 



, J7*2 



'A, 



\. — A miVl A .. TM2fc^2fc _ 

= trace ((QJ) 2k ) 

Since Q is symmetric and J is skew-symmetric, their product QJ is traceless. 
Hence (QJ) 2 = - det(QJ) id. Write u? := det(QJ) = det(Q). Then: 

A Afc =(-\fu? k z»Q^z v = {-\) k u? k 2A 
A Afc =(-l) fc w 2fe trace(id) = {-\) k u? k 2 

As for the coefficients c&, k and CA k , Factum Appendix[UJgives us the relation 
c A fc = ~ 2fcCA fe _i- 

We now plug all this values in l|4.12|l to obtain, for a general function / and 
a quadratic symbol A: 



D 



exp 



, fc=0 



ih 



21, 



' CAfc+l1 (2fc + l)! 



Ei ih 



k=0 



2/.- 



D 



2 k 



\CA k 



(2k)\ 



(6.1) 



f(A) 



We are left with a combinatorics problem: the sequence {cA fc }- Fact El in 
Appendix [HI gives us a recurrence formula: 



k-1 

■E 

J=0 



2k 

27 + 1 



CA, C A , 



(6.2) 
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The sequence is alternating in sign, and the first few absolute values are 1, 2, 
16, 272, 7936 . . .This sequence is called the Zag numbers |18| and they appear 
in the McLaurin expansion of the tangent: 



tana; = 



E 

fe=0 



\CA k 



(2fc+l)! 



„2fe+l 



(6.3) 



To prove this, notice that tana; is the only odd solution to the differential 
equation 



y' = i 



v 



Write a generic solution of the form y = f(x) = Ylk=o ' 



(6.4) 

(2k+l)\ x2k+1 ail< ^ SU ^" 

stitutc it into (|6.4|l . Equating coefficients, we conclude that the sequence atk 
satisfies the same recurrence relation as the sequence CA k (Equation l|6. 2(1 ^1. 
Hence au = |caJ- 

The Zag numbers can be written in terms of the Bernoulli numbers B n : 



l c A fc 



\B 2k \ 



2 2k - 1 2k 

When we use equation l|6.3|) in i|6.1[l we obtain a nice, compact expression: 



B = f(A) 



B 



ihioD 



exp 



2A ifwjD 
— — tan — - — 

irujj 2 



AD 



f{A) 



(6.5) 



Remember that D is the derivative operator that applies to /. If we also 
take the function f(y) = e ev , then D acts simply as multiplication by e. In 
particular, when we consider the time evolution operator: 



B = e - itA ' n 
B 



tu> 

sec — exp 
2 



2A tuj 
— — tan — 

inio 2 



(6.6) 



Equation (|6.6() is derived in ^B] derived in a different manner. For the case 
A = I, it is a well known formula. (See, for instance, 



7 Application: Bohr— Sommerfeld rules. 

We explain in this section the application that caused our original interest in 
this problem. 

In UJ, Colin de Verdiere gives an algorithm which computes the Bohr- 
Sommerfeld quantization rules to all orders in h in the one dimensional case 
N — 1. His method is inspried by Voros [T!5] and a similar method had been 
previously used by Argyres 
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Let H be an operator with symbol H £ C°°(T*R). Bohr-Sommerfeld quan- 
tization rules provide a way to asymptotically compute the spectrum of H. 
Assume H has a regular minimum at a point. Under certain extra assumptions 
on the symbol H (see the eigenvalues of H are given by the solutions E to 

oo 

2nnh = S{E) = ^WS 3 {E) (7.1) 

3=0 

with n G Z. To solve the previous equation, write E = h k Ek as a power 
series in h and substitute it into i|7.1|) to obtain recursive expressions for each 
Ek- This requires knowing the form of S(E), called the semiclassical action, for 
the hamiltonian H. 

Let us consider for simplicity the case where the symbol H does not depend 
on h. It is known that at lowest orders in h: 



S (E)= [ pdx, S 1 (E) = n 

J IE 



'7e 

r-i 



7b denotes the level set H (E) around the minimum of H and {x,p) are the 
natural coordinates in T*K. 

The main result in joj is the following: 

where £ is the parametrization of 7^ by the time evolution 

= Hpdt , dp = ~H x dt 

and P,'^ are universal polynomials in the derivatives of H defined by the symbol 
R a of the resolvent R a = (a — A)^ 1 : 

1 , 



= 77 + > > J 



H ^ ^ (a-HY 

3=1 1=2 K 1 

Looking back at our formula for the symbol for the function of an operator 
l|l.lfl and using it for the function f(y) = (a — y)^ 1 , we see that Lj = 3j/2 + 1 
and the polynomials Pj.i can actually be defined in terms of graphs: 

reduced graphs T 
with I — 1 vertices 
and j ct 



which gives us the following equation for the eigenvalues E of H: 
27r(n-l/2)fi = S Q (J5) 



reduced graphs T 
with B r >0 



7-B 
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Here E is an eigenvalue of H, whereas Et is the number of edges of a graph T. 

Sj(E) = for j > 1 and odd. 82(E) is given by the contribution of 2 
graphs, and 84(E) is given by the contribution of 15 graphs (see AppendixQ. 
However, there is a trick using Stokes' theorem that allow us to express the 
contribution of certain graphs in this expression in terms of others. (This trick 
is used in 4 and in (i for S2, although without the diagrammatic notation.) 
As a consequence, 82(E) can be written in terms of 1 graph and 84(E) can be 
written in terms of 5 graphs (those where every vertex has at least two edges) : 



2%(n-l/2)h = S Q (E) 

~ T 



4= I (H =£ II ) (II 



2! • 6 dE 



~E 



16 



1 d 
2! • 120 dE 



[ (H=^H)dt + -1 

J~1K 



H) 2 dt 



1 



3! -15 \ dE 



IE 



H dt 
H 




H dt 



H) dt 



o(tf) 



(7.3) 



All the integrands in the previous expressions are long polynomials in the 
derivatives of H . The expression would be hard to obtain without the dia- 
grammatic notation. Given a concrete hamiltonian H we could easily program 
a computer to write all the terms in l|7.2|l for that specific operator at higher 
orders in h. 

For the case of a hamiltonian of the form kinetic plus potential energy 
H(x,p) — p 2 /2m + V(x) the contribution of many graphs vanishes, and (|7.3fl 
becomes: 



2„(„-l/2)« = S4E) -£±±f v 



2 7 3 2 



d 

dE 



[V"(x)fdt 



'(x)dt 



d 

dE 



V i4) (x)dt 



O 



The disadvantage of this method is that it does not generalize to the mul- 
tidimensional case N > 1. Cargo et al jj approached this problem in a totally 
different way to obtain a result valid in all dimensions. The symbol of a function 
of an operator plays a role in their derivation, too. 
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A A lemma for calculations with graphs. 



The diagramatic notation introduced in ^ makes equations and derivations 
easier to write. The following lemma is needed in some of those derivations: 

Lemma A.l. Let T be a labeled graph with V vertices and E edges. Let 
D, Ci, . . . , Cy be symbols. Then 



A r (d,...CV) ] k p = ^ \ r >(C u ...,C v ,D) 

r' 

where the sum is taken over all labeled graphs V with V' = V + 1 vertices and 
E' = E + k edges, which are constructed by putting together 

- the labeled graph T ( conserving its labels ), 

- an extra vertex labeled by V + 1, 

- k extra arrows (labeled by E + 1, . . . , E + k) starting from the vertices of 
r and ending at the vertex V + 1 . 

Proof. Write down the definition of both sides in terms of J^ v and check that 
they are equal. □ 

For instance: 



(C^^D)^E 

= C^^D + C—^D^^E 

2 3 

+ C D ^==^ E + C D ^=p* E 

= C^^D E + C ^D^^E 

+ 2 c^^D^^E 



B An alternative derivation of ( 14.41) . 

As we mentioned in the introduction and in [JJJ the existence of a universal for- 
mula like was used by Voros [T§| and Colin de Verdiere 6 (and, indirectly, 
by Argyres [2]). They start by writing a smooth function of an operator in terms 
of the resolvent. Let a £ C and define the resolvent operator R a := (a — A)^ 1 
with symbol R a . Then we use Helffer-Sjostrand's formula ^U] as a spectral 
theorem: 

B = f(A) = -- [ R z chj(z) dxdy 

7T JC 
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Here z = x + iy, / is an almost analytic extension of /, and d-g = d x + id y . This 
allows us to write for the symbol: 

B = -- ( R z dgf(z) dxdy (B.l) 

Hence, finding the symbol of f(A) reduces to finding the symbol of R a ? In 
order to do so, we may write 



k=0 

and, since R a {a — A) = 1 , substitute it into 

R a ★ (a - A) = (a - A) ★ R a = 1 

to obtain recursively the value of each R a ^)- Although simple, this method 
quickly proves intractable. The calculations at order 4 are already too complex 
and we will not find the pattern that l|l.l|) shows. 

But we can also use this idea to prove l|4.4|l in a different way. First, we 
prove it for the function f(y) — {a — y)" 1 : 



Claim. The function 



E\\2 (a- A) v+1 

reduced 
labeled 
graphs r 

satisfies h a {A) * (a — A) = (a — A) * h a {A) = 1. 

The proof is a long combinatorial exercise on calculations with graphs. Therefore 

h a (A) = R a . 

Second, we can substitute (|B.2(I into (|B.1(I to obtain again l|4.4(l . 

C Calculation of cr and St- 

Calculating cr is a combinatorial problem. The following five facts give quick, 
recursive rules for it: 



3 



We could also use Cauchy's integration formula as a spectral theorem: 

B = f(A) = JJ^ f(*)R* 
da 

which will lead to the same results, but it is only valid for analytic functions /. 7 is a path 
around the spectrum of A. 
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Facts: 1. c. = 1 



2. c r = if E is odd. 

This is due to a symmetry property. If we denote by T a labeled graph 
with a numbering of the vertices by 1, 2, • • • V and by i(T) the relabeling 
of the vertices by the permutation (y... 2 i) > then 

Aj( r) = (-l) s A r 

And when we sum over reorderings of the vertices: 

cr = (-l) E cr 

3. If Ti has Vi vertices, T; ^ Tj for i ^ j, and T = T± + ■ ■ -F n is the 
topological sum with V = Vi + • • • V n vertices, then 

Cr= ^!. y! KJ Cri '" Cr " (CJ) 

4. cr does not change if we erase two edges with the same endpoints. 
For instance 

c 1^2=?3^4 == c 1^2 3^4 

5. If p is a vertex in T, denote by T — p the same T with the vertex p and 
every edge starting or ending at p erased. For instance: 

r = • — *-p — - • ■ ■ • 

T—p = • • ? » 

Then: 



( ^ y ^ 1)^ °^ arrows starting at p ^ 

vertices p in T 



Proof. Count the possible reorderings of the vertices by choosing first 
which vertex has label V. □ 

In particular, putting together facts [21 and [21 we get: 

Lemma C.l. If T has a connected component with an odd number of edges, 
then cr = 0. 

These are the reduced, connected graphs with 2 and 4 edges, and for each 
of them the value of cr and Sr- 
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r 


E 


V 


Sr 


cr 


• s- • 


2 


2 


4 


2 


• s- • *- • 


2 


3 


2 


-2 


• g • 


4 


2 


48 


2 


• £ • 5- • 

• s- • s- • 

m s- • 


4 


3 


6 
8 
4 


-2 
6 
2 


i 


> 




4 


4 


8 

2 

4 

2 
4 


8 

8 

-8 









• a- « 


1 

• s- • 


< 

m m 


1 


• s~ • s- • • 


• s- • s- • S- • *- • 


4 


5 


2 
2 

24 


16 

-8 

-24 




< 


1 


< 


• 



D Symbol of a function of an operator at order 
4 in h. 

Using any of the equations derived in S^J we write down the explicit form of all 
the terms of the symbol of a function of an operator up to order 4 in K. The 
data in the table in Appendix IU1 are needed. 
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B = .f(A) 
B = f(A) 

4 

+ 16 



A^^A J" {A) 
2 2! 

A^m A f"( A) 



+ A >■ A ■ 



+ 



24 2! 
A^=^A~* A 



A 



f"(A) 
3! 



A 



A: 



A: 



: A 



+ - {A=*A? + A >A 



+ 4 A *-A *- A ■ 

+ 2 A *A *- A 




f"'(A) 
3! 



f (i) (A) 



4! 



+ 8 A *-A >■ A + A >■ A ■ 



+ 5 ( A j- A ) (A — 

+ 4 A ■* A >■ A ■ 



A A 

A) 



+ 10 (A- 



A) 



2 f {6) (A) 



6! 



0(h 6 ) 



f i5) (A) 
5! 
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